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We report the observation of B~ — > J/ipAp and searches for B~ —+ J/ipT,°p and 5° — > J/ippp 
decays, using a sample of 275 million BB pairs collected with the Belle detector at the T(4S) 
resonance. We observe a signal of 17.2 ± 4.1 events with a significance of 11. la and obtain a 
branching fraction of B(B~ J/ipAp) = 11.6 ± 2.8(stat.)t2;3(sys.) x 10~ 6 . No signal is found for 
either of the two decay modes, B~ — > J/ipYpp and B° — * J/ippp, and upper limits for the branching 
fractions are determined to be B(B~ -> J/y>E°p) < 1.1 x 1CT 5 and B(B° -> J/ippp) < 8.3 x 1CT 7 
at 90% confidence level. 
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We report the observation of the decay mode B~ — > 
J/ipAp, which is a new type of baryonic B decay, B — > 
charmonium + baryons. Evidence for this mode was re- 
ported in previous studies by BaBar 1] and Belle 0, 
with 89 million and 85 million BB pairs, respectively. 

Originally, modes of this type were proposed as a po- 
tential explanation Q for the excess in the low mo- 
mentum region of the inclusive J/ip momentum spec- 
trum for B^J/if + l|| which is not consistent 
with the predictions from nonrelativistic QCD calcula- 
tions 0|. Although the measured branching fraction of 
B~ — > J/ipAp, O(10~ 5 ), is not large enough to explain 
the observed excess, this result stimulated experiments to 
explore new baryonic B decay modes in addition to those 
that had already been observed, namely: B — > charmed 
baryon (+ light anti-baryon and mesons) Q, B — > D 
+ baryons |8J, B — > charmless baryons + light mesons 
0, and B — > baryons + 7 10]. One of the interesting 
features of many three-body baryonic B decays is the 
presence of an enhancement in the threshold region of 
the baryon anti-baryon invariant mass spectrum. Fur- 
ther studies of B — > J /ipAp and related baryonic decay 
modes arc expected to provide additional information on 
the baryon formation mechanism in B decays. 

In this paper, we report the observation of B~ — » 
JjjbAp and searches for B~ — > J /i^Y^p and B° — » J/tpPP 
in a sample of 253 fb _1 containing 275 million BB 
pairs accumulated at the T(AS) resonance with the Belle 
detector a t the KEKB energy asymmetric e + e~ col- 
lider 0. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cerenkov counters (ACC), a barrel- like ar- 



rangement of time-of- flight scintillation counters (TOF) , 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL). These detectors are located inside a su- 
perconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
coil is instrumented to detect Kl mesons and to identify 
muons (KLM). 

Candidates for B~ — > J/ipAp and B° — > Jjvbpp are 
fully reconstructed with all daughters in the final states, 
where we use the decay chains J/i/j — * (I = e,/i) 

and A — > pn~ . Candidates for B" — > J/%pYPp are 
reconstructed partially, where only the daughters from 
the J /if) and A decays plus the anti-proton are recon- 
structed in the final state. In this case, the A origi- 
nates from the S° — ■> A7 decay. Partial reconstruction of 
B~ — > J/ipTrp without the low momentum 7 from the 
E° decay gives better sensitivity than full reconstruction 
because of the low detection efficiency for soft photons. 
As a result, the same selection criteria are applied for 
fully reconstructed B~ — » J/ipAp and partially recon- 
structed B~ -> J/ipE°p. 

Events are required to pass a basic hadronic event 
selection 14] . To suppress continuum background 
(e + e~ — > qq, where q = u,d,s,c), we require the ratio 
of the second to zeroth Fox- Wolfram moments 15] to be 
less than 0.5. 

The selection criteria for the J/ip decaying to are 
identical to those used in our previous papers 0, Il4j . 
J jvb candidates are pairs of oppositely charged tracks 
that originate from within 5 cm of the nominal inter- 
action point (IP) along the beam direction and are pos- 
itively identified as leptons. In order to reduce the ef- 
fect of bremsstrahlung or final state radiation, photons 
detected in the ECL within 0.05 radians of the original 
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e~ or e + direction are included in the calculation of the 
e + e~(7) invariant mass. Because of the radiative low- 
mass tail, the J/tp candidates are required to be within an 
asymmetric invariant mass window: — 150(— 60) MeV/c 2 



GeV/c 2 (-0.32 GeV/c 2 < AM B < 0.20 GeV/c 2 ) and 
Mhc > 5.20 GeV/c 2 are selected for the final analy- 
sis. The signal regions are defined to be 5.27 GeV/c 2 
< M bc < 5.29 GeV/c 2 and |AM B | < 0.03 GcV/c 2 for 



and -0.104 GeV/c 2 < AM B < -0.044 GeV/c 2 for par- 
tially reconstructed B~ -y J/ipYPp, and 5.27 GeV/c 2 
< M bc < 5.29 GeV/c 2 and |AM B | < 0.02 GeV/c 2 for 
B° — > J/tppp, which corresponds to three standard de- 
viations based on the MC simulation. The candidates 
outside of the signal region are used to study background. 

After all selection requirements, about 11% of the 
B~ —y J/ipAp candidates have more than one entry 
per event; this occurs for 20% of the partially recon- 
structed B~ -> J/ipY,°p and 3.8% of the B° -> J/%ppp 
candidates. For these events, the B candidate with the 
best vertex quality is selected. If an event satisfies both 
B~ — > J/tpAp and B° —y J/tppp selection criteria, we 
take it only as a B~ — > J/t/jAp candidate. 

The signal yields are extracted by maximizing the two 
dimensional (2D) extended likelihood function, 
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< M e+e - (7) (M M+M -)-m J/v , <+36(+36)MeV/c 2 , where B~ -► J/^Ap, 5.265 GeV/c 2 < M bc < 5.29 GeV/c 
TOj/y, is the nominal ■//?/> mass 16] . In order to improve 
the momentum resolution of the J /if) signal, a vertex and 
mass constrained fit to the reconstructed J/ip candidates 
is then performed and a loose cut on the vertex fit quality 
is applied. 

In order to identify hadrons, a likelihood Li for each 
hadron type i (i = it, K and p) is formed using in- 
formation from the ACC, TOF and dE/dx. Charged 
tracks that were previously identified as electrons or 
muons are rejected in the hadron identification proce- 
dure. The proton and pion from the A decay are se- 
lected with the requirements of L p /(L p + Lk) > 0.6 and 
L 7r /(L 7r + Lk) > 0.6, which have efficiencies of 94.8% 
and 99.2%, respectively. 

The primary protons from the B~ (B°) decay are se- 
lected with requirements of L p /(L p + Lk) > 0.6 (0.7) 
and L p /(L p + L n ) > 0.9 (0.9). For the primary proton, 
we also apply requirements on dr and dz, the impact pa- 
rameters perpendicular to and along the beam direction 
with respect to the IP, respectively. The efficiency of dr 
and dz requirements is 97.8% (91.2%) for B~ (B°) signal 
modes. 

For A candidates, we impose momentum-dependent re- 
quirements on the dr values for both A daughter tracks, 
the distance between the daughter tracks along the beam 
axis at the A vertex, the difference of azimuthal angles 
between the A momentum and the direction of the A 
vertex from the IP, and the flight length of the A, which 
retain 81.11% of A candidates. The invariant mass of the 
A candidate is required to be within 5 MeV/c 2 (3cr) of 
the A mass. We apply vertex and mass constrained fits 
for the A candidates in order to improve the momentum 
resolution. 

B mesons are reconstructed by combining a J/tp, a 
A(p) and a p candidate. To reduce combinatorial back- 
ground, we impose a requirement on the quality (x 2 ) of 
the vertex fit for the leptons from J/ip and primary pro- 
ton(s) that retains 95.8% (95.1%) o f B~ (B°) si gnal. 

These criteria maximize N S i S / ^/N s i g + N^ s , where 
A^ s ig is the number of expected signal events from sig- 
nal Monte Carlo (MC) samples with assumed branching 
fractions of 1.0 x 10~ 5 for B~ -> J/tpAp and 1.0 x 10~ 6 
for B° — > J/ippp, and A^kg is the number of expected 
background events estimated from the sideband data. 

We identify B candidates using two kinematic variables 
calculated in the center-of-mass system: the beam-energy 



where N is the total number of candidate events, i is 
the identifier of the i-th event, N k and P k are the yield 
and probability density function (PDF) of the compo- 
nent k, which corresponds to the signal and background. 
The signal yields for B~ — > J/ipAp and partially re- 
constructed B~ — y J /ipYPp are simultaneously extracted 
from the B~ — > J/tpAp candidate sample. 

The signal PDFs for all three decay modes are mod- 
eled using a sum of two Gaussians for Mb c and a sum of 
two Gaussians plus a bifurcated Gaussian that describes 
the tail of the distribution for AM B - The PDF param- 
eters are initially determined using signal MC; the pri- 
mary Gaussian parameters are then corrected using the 
B+ -> J/i)K* + {K*+ -> K s ir+) control sample. The 
parameters are kept fixed in the fits to the data. 

The dominant background comes from random com- 
binations of J/tp, A (p), and p candidates. A threshold 
function is used for the M hc PDF. For AM B , we 
take the effect of the kinematic boundary into account 
by using the function 



Pthr{x;x Cl p, c) 



(x - x c ) p e-< x - x J (x > x c ) 
(x < x c ) 



constrained mass (M^ c = \/E^ 



P B ) and the mass 



difference (AM B = M B - m B ) |17j, where E h cam is the 
beam energy, P B and M B are the reconstructed momen- 
tum and mass of the B candidate, and ms is the nom- 
inal B mass. B~ (B°) candidates within \AM B \ < 0.20 



with x — AM B and x c — — (m B — rrij/^p — m\(m p ) — 
m p ), where m B , rrij/^p, m\ (m p ) and m p are the nominal 
masses of the B, J/ip, A (p) and p, respectively. 

For B° —y J/tppp, a cross-feed from B~ — » J/ipAp, 
where the it from the A is undetected, forms a peak in 
the negative AM B sideband region. The PDF for this 
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TABLE I: Summary of the results. Y and b are the fitted signal and total background yields in the signal region, no is the 
observed number of candidate events in the signal region, e (error includes systematic error) is the detection efficiency, Y90 is 
the upper limit on the signal yield at 90% confidence level and B is the branching fraction. 





mode 


Y b 


no 


<%) 


Y 90 


B 


Br 
B~ 
B° 


— > J/ipAp 
-> J/tpTPp 
-> J/ippp 


17.2 ±4.1 0.41 ± 0.09(stat.) 
-1.1 ± 1.7 0.31 ± 0.04(stat.) ± 0.03(sys.) 
-6.1 ± 2.2 0.94 ± 0.10(stat.)t^ ?6( s y s -) 
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FIG. 1: (a) The (M hc , AM B ) scatterplot of B -» J/tpAp candidates and its projections onto (b) AM B with 5.265 GeV/c 2 
< Mbc < 5.29 GeV/c 2 and (c) Mbc with |AMs| < 0.03 GeV/c 2 . The dashed lines and solid boxes indicate the signal regions. 
The curves are the result of the fit as described in the text. 



cross-feed is modeled by a smoothed histogram from the 
B~ -> J/ipAp MC sample. 

In the fit, the value of Nk and the parameters for com- 
binatoric background are allowed to float. Figures Q]and 
El show the (Mbc, AAfg) scatterplots and their projec- 
tions for candidates after all selections are applied. The 
fit results are superimposed on the projections. There 
are sixteen candidate events in the signal region for 
B~ — > J/ipAp, one for B~ — > J/tpYPp and three for 

b° -> j/ippp. 

Table [I] summarizes the maximum-likelihood fit results 
for the signal (Y) and signal-region background (b) yields 
and their statistical errors. For the B~ — > J/ipAp decay 
the fit gives 17.2 ±4.1 signal events with a statistical 
significance of 11. lc. The statistical significance is de- 
fined as a/— 2 ln(£o / £ max ) > where £ max and Co denote 
the maximum likelihood with the fitted signal yield and 
with the yield fixed at zero, respectively. No significant 
signal is found for the B~ — > J/ij)I!Pp and B° — » J/ippp 
decay modes, while the number of cross-feed events in the 
AMb sideband of B° — > J/ippp (9.0 ± 3.2) is consistent 
with the expectation from the observed B~ — > J/ipAp 
signal yield (7.0 ± 1.7). For these two modes, we obtain 
upper limits on the yield at 90% confidence level (I90) 
using the Feldman-Cousins method 0], which takes into 
account the systematic errors due to the uncertainties in 
the signal detection efficiency and the background yield 

m. 



The branching fraction is determined with iV<j/[e x 
N BB x B(J/ip -> 1+1-) x B{A -> pit-)] for B~ -> J/^Ap 
and the partially reconstructed £>~ — > J/i[)Y,°p. For 
5° -> J/^pp we use N s /[e x 7V B5 x B(J/^> -> Z+r)]. 
Here JVs is the signal yield, -ZV B g is the number of -B-B 
pairs. We use the world averages [l(j for the branching 
fractions of B(J/ip -> Z + Z") and fi(A -> jot"). The ef- 
ficiencies (e) are determined from the signal MC sample 
with the same selection as used in the data. A three- 
body phase space model is employed for all three decay 
modes. The fractions of neutral and charged B mesons 
produced in T(4S') decays are assumed to be equal. 

Figure|21demonstrates the consistency between the ob- 
served invariant mass distributions M(J/ijjA), M(J/ipp), 
and M(Ap) of the sixteen B~ —* J/ipAp candidates and 
the phase space distributions obtained from signal MC. 

The sources and sizes of systematic uncertainties are 
summarized in Table [H] The systematic uncertainty on 
the yield is examined by varying each fixed shape param- 
eter by ±1<T. A possible bias in the fitting is studied using 
MC samples; no significant bias is found. The systematic 
uncertainty assigned to the yield is 2.5%. 

The dominant sources of systematic errors on the effi- 
ciency are uncertainties in the three-body decay model, 
A reconstruction, tracking efficiency and particle iden- 
tification. To estimate the error due to uncertainty in 
decay modeling for B~ — > J/ipAp, we subdivide phase 
space into a few bins and recompute the branching frac- 
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FIG. 2: (a) (M bc , AM B ) scatterplot of B° — > J/ippp candidates and its projections on to (b) AM B with 5.27 GeV/c 2 < Mbc < 
5.29 GeV/c 2 and (c) Mb c with |AMb| < 0.02 GeV/c 2 . The dashed lines and the solid box indicate the signal regions. The 
curves are the result of the fit. 
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FIG. 3: Invariant mass distributions for (a) M(J/tpA), (b) M(J/ipp) and (c) M(Ap) for the sixteen signal-region B — * J/ipAp 
candidates. The histograms are phase space distributions from the signal MC sample normalized to sixteen events. 



tion using MC-determined bin- by- bin efficiencies. The 
maximum difference between the nominal value of the 
branching fraction and the values obtained with differ- 
ent choices of bins is assigned as a systematic error. 
For B~ — > J/^Y. p (B° — ► J/ippp), the distribution 
in phase space is unknown and we conservatively assign 
the maximum variation of efficiency among the slices of 
M(JMA), M(JMp) and M(A,p) [M(J/^,p), M(J/^,p) 
and M(p,p)] as a systematic uncertainty. The A recon- 
struction error is determined by comparing the transverse 
proper flight distance distributions for data and MC sim- 
ulation. The discrepancies weighted by distributions of 
the A momentum and the decay length are assigned as 
systematic errors. The uncertainties in the tracking effi- 
ciency are estimated by linearly adding the momentum- 
dependent single track systematic errors. We assign un- 
certainties of 3% per proton, 2% per pion, and 2% per 
lcpton for the particle and lepton identification. 

The systematic errors for the background yield are 
evaluated by varying each of the PDF parameters by its 
statistical error from the fit and by fixing the signal yield 
at zero in the case of B° — > J/ippp to accommodate the 
possibility that an overestimate of the background might 



be causing the large negative signal yield (—6.1 ± 2.2). 

In summary, we have observed a 17.2 ±4.1 event signal 
for B~ — > J/ipAp, with a statistical significance of ll.ler. 
The measured branching fraction is B(B~ — » J/ipAp) 
= 11.6 ± 2.8(stat.)±^;|(sys.) x 1CT 6 . This establishes 
a new type of baryonic B decay, B — » charmonium + 
baryon anti-baryon. The Ap distribution for this decay 
is consistent with a phase space, in contrast to B — > Ap7r 
|9| and other baryonic B decays. No significant signals 
are found for the B~ — > J/ipYfp and B° — > J/ippp 
decay modes. We obtain upper limits on the branch- 
ing fractions of B(B~ — > J/tpT,°p) < 1.1 x 10~ 5 and 
B(B° -> J/ippp) < 8.3 x 1(T 7 at 90% confidence level. 

We thank the KEKB group for the excellent opera- 
tion of the accelerator, the KEK cryogenics group for the 
efficient operation of the solenoid, and the KEK com- 
puter group and the Nil for valuable computing and 
Super-SINET network support. We acknowledge sup- 
port from MEXT and JSPS (Japan); ARC and DEST 
(Australia); NSFC (contract No. 10175071, China); DST 
(India); the BK21 program of MOEHRD, and the CHEP 
SRC and BR (grant No. R01-2005-000-10089-0) pro- 
grams of KOSEF (Korea); KBN (contract No. 2P03B 



6 



TABLE II: Summary of systematic uncertainties (%) in yield and detection efficiency. 



Source 


B~ -> J/ipAp 


B~ -> J/V>S u p 


B u -> Jlihpp 

i i rr 


Uncertainty in yield 


±2.5 






Tracking Error 


±8.6 


±9.8 


±4.8 


PID (proton and pion) 


±8.0 


±8.0 


±6.0 


Lepton ID 


±4.0 


±4.0 


±4.0 


A Reconstruction 


±7.9 


±9.8 




A BR 


±0.8 


±0.8 




MC Statistics 


±0.2 


±0.2 


±0.5 


3-body decay model 


±4.7/-f2.6 


±38.0/-29.1 


+24.1/-18.8 


Total 


+15.7/-19.6 


±41.4/^33.4 


±25.6/-20.7 



01324, Poland); MIST (Russia); MHEST (Slovenia); 
SNSF (Switzerland); NSC and MOE (Taiwan); and DOE 
(USA). 



* Present adderess: Physics Department of Sichuan Uni- 
versity, Sichuan Province, China 
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